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ABSTRACT 

In a search with the Parkes radio telescope of 56 unidentified Fermi-L AT gamma-ray sources, we have 
detected 11 millisecond pulsars (MSPs), 10 of them discoveries, of which five were reported in iKerr et alj 
(120121) . We did not detect radio pulsations from another six pulsars now known in these sources. We describe 
the completed survey, which included multiple observations of many targets done to minimize the impact of 
interstellar scintillation, acceleration effects in binary systems, and eclipses. We consider that 23 of the 39 
remaining sources may still be viable pulsar candidates. We present timing solutions and polarimetry for 
five of the MSPs, and gamma-ray pulsations for PSR J1903-7051 (pulsations for five others were reported 
in the second Fermi-L AT catalog of gamma-ray pulsars). Two of the new MSPs are isolated and five are 
in > 1 d circular orbits with O.2-O.3M 0 presumed white dwarf companions. PSR J0955-6150, in a 24d 
orbit with a « 0.25 M 0 companion but eccentricity of 0.11, belongs to a recently identified class of eccentric 
MSPs. PSR J1036—8317 is in an 8 hr binary with a > O.14M 0 companion that is probably a white dwarf. 

PSR J1946-5403 is in a 3 hr orbit with a > 0.02 M 0 companion with no evidence of radio eclipses. 

Subject headings: gamma-rays: stars — pulsars: individual (PSR J0955-6150, PSR J1012-4235, 

PSR J1036—8317, PSR J1903-7051, PSR J1946-5403) 


1. INTRODUCTION 

The Large Area Telescope (LAT; lAtwood et al.l 12009! ) on 
the Fermi Gamma-ray Space Telescope is a superb instrument 
with which to study rotation-powered pulsars. Since mid- 
2008, it has been used to detect more than 170 pulsars at en¬ 
ergies above 0.1 GeV (e.g.. lAbdo et al 11201 Obi pi Fermi has 
identified millisecond puls ars ( MSPs ) as a ubiquitous class of 
gamma-ray sources (e.g., lAbdo et aljl2009al) . These old neu¬ 
tron stars, spun up by accretion from an evolved companion 
(lAlpar et all 19821) . are a relatively local and isotropically dis¬ 
tributed population, that make up only 10% of the identified 
population of pulsars in the Galactic disk. However, about 
half of known gamma-ray pulsars are MSPs. 

Nearly half of the known gamma-ray MSPs were discov¬ 
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ered as radio objects in undirected (“all sky”) surveys prior 
to the launch of Fermi. Gamma-ray pulsations were sub¬ 
sequently detected with the aid of rotational ephemerides 
obtained from radio timing observations. Many slowly- 
rotating pulsars have been discovered via direct periodic- 
ity searches of sparse gam ma-ray photons (e.g., Ab do et all 
2009b: IPletsch et al.l2012hl) . but so far this has not been possi¬ 
ble unbiasedly for binary MSPs. However, the LAT has led to 
the discovery of MSPs in a di fferent, and prolific fashion. In 
the three LAT source cat alogs (lAbdo et al 1201 Oat iNolan et alJ 
1201 2t lAcero et alJ 120151) . there are hundreds of unidentified 
sources, many of which have spectral characteristics typical 
of pulsars. Radio searches of many of these have turned up 
dozens of MSPs so far, and once the radio ephemerides have 
been obtain ed, gamma-ray pulsa t ions have almost always fol¬ 
lowed (e.g., Co gnard et al.ll201 ltlBarr et al.ll2013l) . 

Using the CSIRO Parkes telescope in 2009, we discovered 
five MSPs in a radio survey of 14 unidentified LAT sources 
(IKerr et al .1120121) . In an extension of that survey we have dis¬ 
covered five more MSPs. Here we present the completed sur¬ 
vey, and report on radio timing, polarimetric, and gamma-ray 
studies of some of the MSPs. 

2. OBSERVATIONS, ANALYSIS, AND RESULTS 

In this Section we describe the radio searches of unidenti¬ 
fied gamma-ra y sou rces that we performed at the Parkes tele¬ 
scope (Section IXT1) . the s ensit ivity of the survey and relevant 
selection effects (Section 12.21) . radio timing and pol ari metri c 
observations of the pulsars discovered (Secti onsl2.3l and l2~4l) . 
gamma-ray results for one of th em (Section [2.5l i. and some 
X-ray observations (Section [2~6l i. 

2.1. Parkes Radio Searches 
2.1.1. Initial Searches of 1FGL Sources 

iKeith et akl (1201 ll) used a digital filterbank at Parkes to 
search 11 unidentified sources from the first Fermi LAT cat- 
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alog (1FGL: lAbdo et all 1201 Oal) . Two MSPs and one slow 
pulsar were discovered in single observations of each target 
at a central frequency of 1.4 GHz. Subsequently, one of these 
MSPs was found to be associated with the corresponding LAT 
source. 

At nearly the same time, in late 2009, we used an analog 
interbank at Parkes to search 14 unidentified 1FGL sources 
dKerr et alj|2012i) . The single observations of these targets re¬ 
sulted in the detection of six MSPs, five of them discoveries. 
However, confirmation of some of these MSPs was not easy: 
the search observations lasted for 1-2 hr each, but some of the 
pulsars were not detected in equivalent initial confirmation at¬ 
tempts, owing to the effects of interstell ar sc intillation. These 
and other selection effects (see Section l±2l for details) led us 
to search some promising unidentified LAT sources repeat¬ 
edly. 

2.1.2. Repeated Searches of Unidentified LAT Sources 

Our subsequent searches used the same equipment and 
methods as lKerr et alJ (120121 ). In brief, total-power measures 
from the central beam of the 20 cm multibeam receiver were 
filtered into 512 contiguous 0.5 MHz-wide channels centered 
on 1390 MHz and sampled 8000 times per second, then digi¬ 
tized with 1-bit precision and written to disk for off-line anal¬ 
ysis. Individual integration times were about 1 hr, and each 
LAT source was observed between one and nine times, de¬ 
pending on the then-perceived quality of the source and tele¬ 
scope availability (Table [I]). 

All the data were analyzed using PRESTO (lRansomll2001l) . 
which implements standard pulsar search techniques includ¬ 
ing radio-frequency interference excision and optimization 
for signals with changing apparent spin periods caused by or¬ 
bital motion. Finite sampling time and smearing from disper¬ 
sive propagation delays within finite-width filterbank chan¬ 
nels both unavoidably degrade sensitivity to pulsed signals. 
Additional smearing can result from use of an incorrect dis¬ 
persion measure (DM) to remove the delays between chan¬ 
nels. We dedispersed each observation using a set of trial 
DMs such that this effect was negligible. This was done up 
to twic e the maximum DM predicted by the ICordes & Laziol 
(120021) model for the corresponding line of sight (see TableQ). 
The maximum acceleration searched for corresponded to sig¬ 
nals drifting by ±200 /iif, bins in the Fourier domain, where 
rih is the largest harmonic at which a signal is detected (up 
to 16 harmonics were summed, in powers of two). This was 
parameterized by zmax = 200 within PRESTO ( Ransom et al.l 
120021 ) . 

In the 2009 searches (Section l2.1.1b we detected six MSPs 
in single observations of 14 LAT sources (TableQ] above the 
horizontal dividing line). In 2010-2011 we re-observed seven 
of the remaining eight 1FGL sources (Table Ol above the di¬ 
viding line), but no new pulsars were detecteco 

The 1FGL catalog was_based on 11 months of LAT data. 
The 2FGL catalog dNolan et al.ll2012l) was based on two years 
of data. In 2012 we selected new search targets, based on 
a three-year source list developed by the LAT collaboration 
but never published (the subsequent 3FGL catalog is based 
on fo ur years of data and a differ ent pipeline: lAcero et al.l 
120151) . As for the 1FGL searches dKerr et all 120121) . we re¬ 
stricted ourselves to non-variable southern sources with no 

14 One of these sources is now known to harbor PSR J1227-4853, which 
recently transitioned to a radio-emitting state tRov et ani2015l and see Ta- 
ble[2>. 


plausible known blazar counterparts and with a LAT posi¬ 
tional uncertainty (95% confidence level error radius) <!', to 
fit within the 1.4 GHz primary beam of the Parkes telescope. 
The remaining sources were classified by visual inspection of 
the gamma-ray spectral energy distribution to pick out candi¬ 
dates with spectral shape resembling those of known pulsars, 
which typically have exponentially cut-off power-law spectra 
(lAbdo et al.l2013l) . Spectral modeling and source localization 
is more difficult for LAT sources near the Galactic plane, and 
we only considered those with \b\ > 4°. 

The new target set consisted of 49 sources, each observed 
between one and seven times, for 122 integrations in the ag¬ 
gregate. Seven of the 49 targets had also been observed in 
their prior 1FGL incarnation (italicized in TableQ). Among 
the remaining 42, we discovered five MSPs (TableQ). 

Overall, we searched 56 individual LAT sources in our 
Parkes survey using the analog filterbank system. We detected 
11 MSPs, of which 10 were discoveries (TableQ). In fact, ra¬ 
dio and/or gamma-ray pulsars are now known in 17 of those 
56 sources (Table [2]). We discuss the sources, in particular 
which might still be good pulsar candidates, in Section Q] 

Three of the five new MSPs were detected unbiasedly 
in all their search observations (i.e., without prior knowl¬ 
edge of their DM or approximate spin period). However, 
PSR J1946-5403 was detected in only three of six observa¬ 
tions, and PSR J1036—8317 was detected in only the third of 
three search observations. We now consider the reasons be¬ 
hind these failures to consistently detect some new pulsars. 

2.2. Sensitivity and Selection Effects 

The nominal sensitivity of our survey to a P = 2 ms pulsar 
with a duty cycle of 25% and DM < 40pc cm -3 , for a typical 
integration time of 1 hr (Table Q), and for a sky temperature 
corresponding to the average at the locations of the targets 
searched, was 0.2 mJy — provided that the dilution of power 
in the Fourier domain caused by orbital motion in a putative 
binary was ideally corrected by our acceleration search analy¬ 
sis. However, two of the newly discovered MSPs have orbital 
periods I), su ch that the discovery integrations were 0.1-0.3 1), 
(Section l2.3l i. for which this idealized correction breaks down 
badly; how badly, for a given system, depends on the observed 
orbital phases (see Figure Q] and lJohnston & Kulkarnili~i~99 ll 
iBagchi et al.ll201 3l) . 

As well, a considerable fraction of LAT-selected MSPs is 
in eclipsing binary systems (e.g.. lHessels et al]l2011l) . Multi¬ 
ple observations of a promising source can thus help combat 
these factors affecting detectability. In addition, the flux den¬ 
sity received from many MSPs is severely affected by inter¬ 
stellar scintillation (see, e.g., Figure 1 of iLevin et al.li2013l) . 
As we now show, this has played a very significant role in 
our survey, and the early realization of the magnitude of some 
of these selection effects is what led us to do multiple ob- 
servatio ns of m any unidentified LAT sources starting in 2010 
(Sec tion !2. 1.2b . 

2.2.1. Interstellar Scintillation and Detection Statistics 

In order to acquire the detections required to deter mine 
the timing solution for PSR J1514-4946 (Section 12.3b . we 
used 100 hr of Parkes telescope time in 79 observations on 
54 days spread over 2 yr. In only 27 of those observations 
did we detect the pulsar in a relatively unbiased manner — by 
dedispersing the raw data using the known DM, performing 
an acceleration search, and looking for a signal with period 
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Table 1 

Radio Searches of FGL Sources at Parkes: Observations 


Name 2 

R.A. b 

Decl. b 

1 

b 

Integration time 

DM max c 


(J2000.0) 

(J2000.0) 

(deg) 

(deg) 

(min) 

(pc cm -3 ) 

3FGL J0101.0-6422 

01 h 00 m 58 s 

-64°24'06" 

301.2 

-52.7 

60 

270 

3FGL J0602.8-4016 

06 h 03 m 05 s 

-40°11'04" 

246.8 

-25.9 

90, 120 

270 

3 FGL J0933.9-6232 

09 h 34 m 00 s 

-62°30'17" 

282.2 

-7.8 

120, 60, 60, 60, 60 

435 

3FGLJ1035.7-6720 

10 h 36 m 10 s 

-67°21'05" 

290.4 

-7.8 

120, 86, 136, 60, 60, 60 

435 

3FGL J1227.9-4854 

12 h 27 m 50 s 

-48°51'57" 

298.9 

13.8 

120, 65, 60 

307 

3FGL J1231.6-5113 

12 h 31 m 49 s 

-51° 18'49" 

299.8 

11.4 

120 

270 

3FGL J1514.2-4947 

15 h 14 m 06 s 

-49°45'33" 

325.2 

6.8 

120 

270 

3FGL J1624.2-4041 

16 h 24 m 07 s 

-40°41'20” 

340.6 

6.2 

120, 120, 120, 60, 60, 60 

435 

3FGL J1658.4-5323 

16 h 58 m 43 s 

-53° 17'43” 

335.0 

-6.6 

83 

270 

3FGLJ1744.1-7619 

17 h 44 m 02 s 

-76°20'25" 

317.1 

- 22.5 

41,90, 80, 72,71,78,60,60,60 

192 

3FGL J1747.6-4037 

17 h 47 m 29 s 

-40° 36'07" 

350.2 

- 6.4 

86 

270 

3FGL J1902.0-5107 

19 h 02 m 05 s 

-51°09'45” 

345.6 

- 22.4 

70 

270 

3FGL J2039.6-5618 

20 h 39 m 30 s 

-56°20'45" 

341.2 

- 37.1 

75, 112, 60 

270 

3FGL J2241.6-5237 d 

22 h 41 m 52 s 

-52°37'38" 

337.4 

- 54.9 

120 

270 

3FGL J0133.0-4413 

01 1 '33 m 27 s 

-44°08'29" 

279.2 

-71.0 

60, 60, 60, 60, 60, 60 

77 

3FGL J0216.1-7016 

02 h 14 m 04 s 

-69°52'06" 

292.9 

-45.6 

60 

115 

3FGL .10744.8-4028 

07 h 44 m 59 s 

-40°29'49" 

254.6 

-8.0 

60 

717 

3FGL J0802.3-5610 

08 h 02 m 46 s 

-56° 15'32" 

270.0 

-13.2 

60, 60, 35, 35, 35 

627 

3FGL J0933.9-6232 

09 h 34 m 02 s 

-62°31'34" 

282.2 

-7.9 

60, 35 

627 

3FGL J0940.6-7609 

09 h 40 m 45 s 

-76°09'34" 

292.2 

-17.4 

60 

269 

3FGL J0940.7-6102 

09 h 40 m 57 s 

-61°05'10" 

281.9 

-6.3 

60 

653 

3FGL J0954.8-3948 

09 h 55 m 03 s 

- 39049 ' 25 " 

269.9 

11.5 

60, 60, 60, 60, 60 

371 

3FGL J0955.6-6148 

09 h 55 m 39 s 

-61°48'36" 

283.7 

-5.7 

60, 60 

653 

3FGL J1012.0-4235 

10 h 12 m 07 s 

-42°35'01" 

274.2 

11.2 

60 

371 

3FGLJ1025.1-6507 

10 h 25 m 00 s 

-65°07'22" 

288.3 

-6.5 

60 

755 

3FGL J1035.7-6720 

10 h 36 m 09 s 

-67°22'17" 

290.4 

-7.8 

45, 60, 60 

563 

3FGL J1036.0-8317 

10 h 36 m 20 s 

-83°17'01" 

298.9 

-21.5 

60, 60, 60 

192 

3FGL J1057.7-6624 

10 h 58 m 43 s 

-66°21'59" 

292.0 

-5.9 

60, 60, 60, 60 

755 

3FGLJ1136.6-6826 

ll h 36 m 47 s 

-68°25'37" 

296.1 

-6.5 

60 

755 

3FGL J1227.9-4854 

12 h 27 m 42 s 

-48°53'28" 

298.9 

13.8 

60, 60 

371 

3FGL J1231.6-5113 

12 h 31 m 34 s 

-51°12'31" 

299.8 

11.5 

60 

435 

3FGL J1238.3-4543 

12 h 38 m 19 s 

- 45043 ' 18 " 

300.5 

17.1 

60 

269 

3FGLJ1306.8-4031 

13 h 06 m 52 s 

-40°32'35" 

306.1 

22.2 

60 

192 

3FGLJ1311.8-3430 

13 h ll m 46 s 

-34°29'19" 

307.7 

28.2 

60, 58, 60, 58 

154 

3FGL J1325.2-5411 

13 h 25 m 15 s 

-54°11'13" 

307.9 

8.4 

60 

627 

3FGL J1326.7-4727 

13 h 26 m 40 s 

-47°27'52" 

309.1 

15.0 

60, 60 

371 

3FGL J1417.5-4402 

14 h 17 m 30 s 

-44°02'40" 

318.9 

16.1 

60 

307 

3FGL J1417.7-5026 

14 h 17 m 39 s 

-50°25'33" 

316.7 

10.1 

60 

627 

3FGL J1518.2-5232 

15 h 18 m 27 s 

-52°33'57" 

324.3 

4.1 

60, 60, 60, 60 

1044 

3FGL J1536.3-4949 

15 h 36 m 29 s 

-49°49'45" 

328.2 

4.8 

37, 37, 16, 60, 60, 60 

883 

3FGL J1539.2-3324 

15 h 39 m 15 s 

-33°25'42" 

338.7 

17.5 

60, 15, 52, 60 

269 

3FGL J1603.7-6011 

16 h 03 m 44 s 

-60°11'12" 

324.8 

-5.7 

60 

883 

3FGL J1617.4-5846 

16 h 17 m 28 s 

-58°46'19" 

327.0 

-5.9 

60 

883 

3FGL J1624.2-4041 

16 h 24 m 09 s 

-40°40'23" 

340.6 

6.2 

60, 60 

883 

3FGL J1702.8-5656 

17 h 02 m 33 s 

-56°54'54" 

332.4 

-9.2 

60, 60, 60, 60 

627 

3FGL J1717.4-5157 

17 h 17 m 35 s 

-51°57'58" 

337.7 

-8.1 

60 

755 

3FGL J1736.2-4444 

17 h 36 m 13 s 

-44°44'50" 

345.5 

-6.7 

60 

883 

3FGLJ1744.1-7619 

17 h 44 m 11 s 

-76°20'29" 

317.1 

-22.5 

60 

192 

3FGL J1753.6-4447 

17 h 53 m 38 s 

-44°47'50" 

347.1 

-9.4 

60 

627 

3FGLJ1803.3-6706 

18 h 03 m 25 s 

-67°07' 14" 

326.9 

-20.4 

60 

230 

3FGL J1808.3-3357 

18 h 08 m 22 s 

-33°56'05" 

358.1 

-6.7 

56, 60, 60, 60 

755 

3FGL J1831.6-6503 

18 h 31 m 45 s 

-65° 05'19" 

329.9 

-22.4 

60 

192 

3FGL J1903.6-7052 

19 h 02 m 43 s 

-70°53'49" 

324.3 

-26.4 

60 

269 

3FGL J1946.4-5403 

19 h 46 m 24 s 

-54°02'46" 

343.9 

-29.6 

60, 60, 60, 60. 60. 54 

154 

3FGL J1959.8-4725 

19 h 59 m 57 s 

-47°26'32" 

351.8 

-30.9 

58, 60, 60, 35, 60 

307 

3FGL J2039.6-5618 

20 h 39 m 51 s 

-56°20'26" 

341.2 

-37.2 

60, 60, 60, 60 

115 

3FGL J2043.8-4801 

20 h 43 m 49 s 

-48°00'45" 

351.7 

-38.3 

60, 60 

115 

3FGLJ2112.5-3044 

21 h 12 m 36 s 

-30°42'37" 

14.9 

-42.4 

56, 60, 60, 60, 60, 60 

103 

3FGLJ2131.1-6625 

21 h 31 m 06 s 

-66°24'42" 

326.7 

-40.3 

60, 60 

115 

3FGLJ2133.0-6433 

21 h 33 m 30 s 

-64°31'58" 

328.7 

-41.3 

60, 60, 60, 60, 60, 60, 60 

103 

3FGL J2200.0-6930 

22 h 00 m 00 s 

-69°31'37" 

321.3 

-40.9 

60 

115 

3FGL J2220.6-6833 

22 h 20 m 24 s 

-68°32'22" 

320.8 

-43.0 

60, 60, 60 

115 

3FGL J2333.0-5525 

23 h 33 m 04 s 

-55°25'32" 

324.2 

-58.3 

60, 60 

87 


Note. — Boldfaced entri es denote observations with detection of MSPs. Discoveries in single-observation searches of the first 14 entries (listed above the horizontal line) were 
reported in lKerr et alj 12013) . Italicized entries denote sources re-observed in 2012 at the improved locations listed below the horizontal line. 

a The names given are of the 3FGL sources closest to our pointing locations (the offset between the pointing positions given here and the 3FGL positions are provided in Table[2l. 
k Parkes telescope pointing position. 

c Maximum trial dispersion measure used in our analysis of the respective data set(s), corresponding approximately to twice the maximum DM predicted by the lCordes & Lazio] 
nm model for the corresponding line of sight. The first observation of each source above the horizontal line was analyzed with DM max = 270 pc cm' 3 (Kerr et al|2012l) . 
d MSP discovered independently bv lKeith et alj 1201 lb . 
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Table 2 

Radio Searches of FGL Sources at Parkes: Source Information 


3FGL name 3 

a e b 

(deg) 

r95 c 

(deg) 

Class' 1 

Sig e 

(o') 

Curve f 

O) 

Var 8 

Spectrum 

notes* 1 


J0101.0-6422i 

0.03 

0.04 

PSR 

26.8 

6.8 

55 

i 

CpR 

i 

J0133.0 4413 

0.12 

0.09 

bll 

7.7 

0.5 

48 

4 

lhr 

6 

J0216.1-7016 

0.44 

0.12 

bcu 

6.0 

0.5 

49 

4 

Id 

1 

J0602.8 4016 

0.10 

0.04 

bcu 

15.3 

1.8 

49 

4 

lihr 

2 

J0711.S '1028 

0.04 

0.08 

bcu 

7.4 

1.9 

39 

3 

Id 

1 

J0S02.3 5610 

0.11 

0.10 


12.3 

3.6 

44 

4 

Id 

5 

J0933.9-6232 

0.01 

0.04 


22.0 

8.8 

59 

1 

CPr 

2 

J0940.6—7609 

0.01 

0.10 


10.0 

1.1 

48 

3 

Id 

1 

J0910.7 6102 

0.05 

0.21 

bcu 

6.4 

2.7 

51 

3 

pr 

1 

J0954.8-3948 

0.04 

0.09 


18.8 

3.9 

51 

3 

D* 

5 

J0955.6-6148 k 

0.00 

0.11 

psr 

7.0 

2.2 

39 

2 

cr 

2 

J1012.0-4235 k 

0.02 

0.07 

psr 

8.8 

3.6 

38 

3 

?P r 

1 

J1025.1-6507 

0.02 

0.09 


7.3 

2.7 

52 

3 

Id 

1 

J1035.7-6720 1 

0.06 

0.04 


30.4 

8.3 

47 

1 

CpR 

3 

J1036.0-8317 k 

0.01 

0.12 

psr 

6.6 

2.3 

42 

3 

?P 

3 

J1057.7 6621 

0.10 

0.06 


9.1 

1.2 

29 

4 

LD 

4 

J1136.6 6826 

0.01 

0.11 

bcu 

7.7 

0.7 

43 

4 

LD 

1 

J1227.9-4854" 1 

0.04 

0.04 

psr 

38.6 

4.6 

74 

3 

Dh* 

2 

J1231.6—5113 

0.01 

0.11 

13.4 

5.6 

46 

1 

CPr 

1 

J1238.3 4543 

0.01 

0.08 

bcu 

8.0 

0.3 

61 

4 

?lh 

1 

J1306.8-1031 

0.01 

0.07 


14.9 

0.3 

43 

4 

ID 

1 

J1311.8-3430" 

0.02 

0.02 

PSR 

62.3 

9.0 

52 

1 

CR 

4 

J1325.2—5411 

0.00 

0.10 


7.0 

2.5 

31 

3 

?ld 

1 

J 1326.7-4727 

0.02 

0.05 

glc 

11.3 

6.1 

66 

1 

Cpr 

2 

J1417.5—4402 

0.00 

0.06 


12.5 

1.4 

54 

2 

cd 

1 

J1117.7 5026 

0.02 

0.11 


6.1 

0.4 

40 

4 

Ld 

1 

J1514.2-4947J 

0.03 

0.02 

PSR 

39.4 

9.3 

35 

1 

CpR 

1 

J1518.2—5232 

0.04 

0.07 


10.3 

3.6 

62 

2 

cr 

4 

J1536.3-4949 0 

0.02 

0.02 

psr 

60.0 

8.5 

51 

3 

pRh 

6 

J1539.2-3324 

0.02 

0.04 


19.2 

9.6 

57 

2 

cPr 

4 

J 1603.7 6011 

0.01 

0.09 


5.0 

3.6 

50 

4 

?h 

1 

J1617.1 5816 

0.00 

0.08 

fsrq 

15.8 

1.6 

128 

5 

LDV 

1 

J1624.2-4041 1 

0.02 

0.04 


19.2 

7.3 

50 

1 

cpR 

2 

J1658.4-5323) 

0.10 

0.06 

PSR 

16.9 

6.5 

38 

1 

Cr 

1 

J1702.8-5656 

0.05 

0.04 


28.8 

6.1 

58 

3 

rdh* 

4 

J1717.1 5157 

0.02 

0.09 

fsrq 

12.2 

2.5 

383 

5 

LDV 

1 

J 1736.2-4444 

0.00 

0.05 

glc 

16.9 

3.3 

46 

2 

cd 

1 

J1744.1-7619 1 

0.01 

0.03 


32.8 

9.9 

51 

1 

CPR 

1 

J1747.6-4037J 

0.05 

0.07 

PSR 

10.0 

2.8 

46 

2 

7c 

1 

J1753.6-4447 

0.01 

0.08 


11.1 

4.0 

41 

1 

cpr 

1 

J1803.3-6706 

0.01 

0.08 


10.4 

1.8 

43 

3 

ldh 

1 

J1808.3-3357 

0.03 

0.09 


8.7 

4.4 

50 

1 

cpr 

4 

J1831.6-6503 

0.02 

0.10 


8.5 

4.5 

35 

2 

cPr 

1 

J1902.0-5107J 

0.04 

0.04 

PSR 

28.9 

6.2 

50 

1 

CR 

1 

J1903.6-7052 k 

0.08 

0.05 

PSR 

16.3 

1.9 

52 

2 

cd 

1 

J1946.4-5403 k 

0.01 

0.06 


20.2 

6.8 

39 

1 

Pi' 

6 

J1959.8 4725 

0.02 

0.03 

bcu 

19.2 

3.9 

49 

4 

pihr 

5 

J2039.6-5618 

0.04 

0.04 


25.3 

5.1 

34 

1 

CpR 

4 

J2043.8—4801 

0.01 

0.08 


9.5 

3.0 

35 

2 

cr 

2 

J2112.5-3044 

0.02 

0.04 


30.1 

7.6 

51 

1 

CpR 

6 

J2131.1—6625 

0.01 

0.11 


10.1 

3.1 

52 

2 

rl 

2 

J2133.0-6433 

0.05 

0.10 


10.1 

4.7 

49 

2 

Pr 

7 

J2200.0—6930 

0.02 

0.13 


9.5 

0.6 

41 

3 

ld 

1 

J2220.6 6833 

0.02 

0.11 


5.7 

1.1 

50 

4 

lh 

3 

J2241.6-5237P 

0.03 

0.03 

PSR 

51.7 

12.6 

60 

1 

CpR 

1 

J2333.0-5525 

0.00 

0.08 


10.7 

2.5 

42 

2 

cr 

2 

a Boldfaced names denote 17 3FGL sources with now known associated radio and/or gamma-ray pulsars; struck-through 

names indicate 16 

i sources 

that we believe are 

no longer 


viable pulsar candidates (see next-to-last column). 

b Offset between 3FGL position and pointing position in Table [T| (where there are two of the latter, only the second is listed here); the Parkes beam half-width at half-maximum is 

0 . 12 °. __ 

c Size of 3FGL source error box (95% confidence level semi-major axis; all 3FGL properties listed here are taken from the 3FGL catalog. [Acero et all2015l) . 

^ Classification from 3FGL pipeline. “PSR” is a pulsar with LAT pulsations; “psr” is a positionally coincident pulsar so far without LAT pulsations; “bll” is a BL Lac; “bcu” is an 
unclassified blazar; “glc” is a globular cluster; “fsrq” is a flat-spectrum radio quasar. 
e 3FGL source significance. 

^ Significance of curvature of 3FGL source spectrum when fit to a log-parabolic model. 

8 Variability index of source (> 73 indicates variability at the > 99% C.L.). 

^ See Section [331 for a description of these classifications and characteristics. 

1 Number of observations of each target at position closest to 3FGL source (from TablelTl. A 0 offset from it. 

J Radio MSP discoveries from our Parkes survey, first reported in lKerr et aljf2012l) . 
k First reported in this work. 

* Pulsar discovered via gamma-ray pulsations (H. J. Pletsch, private communication). 
m Intermittently radio-emitting MSP t|Roy et al|2015l) . 

n MSP in 93-minute orbit discovered via gamma-ray pulsations fPletsch et alj20123) , subsequently detected at the GBT iRay et al|2013l) . 

° Radio MSP discovered at the GMRT i Ray et al. 20T%. 

P Radio and gamma-ray MSP, discovered at Parkes 1 Keith et al|2011|) . 
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Figure 1. Pulse profiles (repeated twice) as a function of time during the 
1 hr observation times, and integrated (at top), for each of the three unbiased 
detections of PSR J1946—5403, folded according to the best values of spin 
period and period derivative determined by the search software. In the other 
three search observations of this source (Table[I} the pulsar was not detected 
unbiasedly, because of the dilution of power in the Fourier domain caused by 
the rapid non-linear change in projected orbital velocity for this 3 hr binary 
system. A constant-acceleration search such as the ones we performed cannot 
correct for this effect well enough. 


3.589 ms. In several of those observations the pulsar was not 
detectable without prior knowledge of the DM and approxi¬ 
mate period. A pulsar like PSR J1514-4946 is therefore dis¬ 
coverable in a search like ours at Parkes less than 1/3 of the 
time! This is due to its small intrinsic flux density combined 
with very large modulation owing to propagation through the 
dynamic and inhomogeneous interstellar medium (see Fig¬ 
ure [2}. 

This large modulation of observed flux density is pri¬ 
marily caused by strong diffractive interstellar scintillation. 
This causes the detected pulsar signal strength in the time- 
frequency plane to form patches (or “scintles”) of charac¬ 
teristic size A v and At. For a given observing frequency, 
A// increases strongly with decreasing pulsar distance d. 
while At increases particularly with decreasing pulsar veloc- 


itv in the nlane of the skv V i (see, e. 2 .. 

Johnston et al. 1998!; 

Nicastro et al. 2001; Lorimer & Kramer 

2005)“ When Av or 


Ar are large compared to the observing bandwidth and inte¬ 
gration time, deep fluctuations in observed flux density result. 
While we have not measured Av or At for PSR J1514-4946, 
the observations (e.g.. Figure 0 are consistent with large val¬ 
ues for both (e.g., > 100 MHz and > 1 hr, respectively), ac¬ 
counting for the very large observed modulations. 

The measured flux density of PSR J1658-5324 is affected 
by scintillation even more, ranging over a factor > 50. It was 
detected in a relatively unbiased manner only 70% of the time, 
and less than 1/2 of the time in the absence of prior DM and 
period information. PSR J0101-6422, another MSP discov¬ 
ered in the first phase of our survey, for which we have good 
statistics, is detectable 80% of the time dKerr et all 120121) . 
These detection statistics are based on t ypical 1 hr individual 
Parkes timing observations (Section l2.3l ). 

For PSR 11903-7051, the observed flux densities range 
over a factor > 10. PSR J1902-5105 varies less: its recorded 
flux density has a standard deviation of 40% of the mean. 


These MSPs are reasonably bright, with relatively narrow 
pulse profiles, and were always detected in typical 20 min 
Parkes timing observations. For PSR J1747-4036, with a 
large DM = 153 pc cm -3 , the effects of scintillation modulate 
the flux density that we record by only 25% about the mean, 
and it w as de tected 100% of the time at Parkes and the GBT 
(Section l2.3l ). 

Away from the Galactic plane, putative pulsar counterparts 
to unidentified gamma-ray sources are largely expected to be 
relatively nearby MSPs, which as a c lass have relatively small 
space velocities (e.g.. iGonzalez et alTI 1201 11) , Depending on 
observing parameters (frequency, bandwidth, and integration 
time), these characteristics can make them particularly sus¬ 
ceptible to deep flux density fluctuations, in turn with im¬ 
portant implications for radio searches. Scintillation modu¬ 
lations in the strong regime have exponential statistics (e.g., 
iRicketti fl 990t) . with median flux density measurements less 
than the mean. Two of the six MSPs fo r wh ich we have 
already obtained timing solutions (Section 12.31 ) display flux 
densities that vary by a factor of > 40, with small median 
values and, given scintillation statistics, not detectable in an 
unbiased manner most of the time for the parameters of our 
Parkes survey. 

The foregoing strongly suggests that an unidentified LAT 
source that is a good pulsar candidate should be searched re¬ 
peatedly before much of a statement can be made about the 
likelihood of it being a radio MSP beaming towards the Earth. 
Eight of the 10 MSPs that we discovered were detected on 
their first observations (TableQ]), but as the preceding account 
illustrates, in some cases (particularly for PSRs J1514-4946 
and 11658-5324) this was fortuitous. Based on our empirical 
evidence, we judge that every promising LAT source should 
be observed a minimum of three or four times in a survey such 
as the one we did at Parkes before it can reasonably be con¬ 
sidered searched in an average sense. 

As indicated in Table [2] we consider that a number of our 
survey targets remain promising pulsar candidates, and at 
least some shou ld be re -sear ched in radio. We discuss this 
more in Sections [ 33 ] and [U6] 

2.3. Radio Timing 

We began timing observations of all 10 MSPs immediately 
following their discovery. In Table [ 3 ] we present initial param¬ 
eters for the four most recent discoveries. We have determined 
phase-connected rotational ephemerides for the six remaining 
MSPs. That for PSR J0101-6422 was reported in I Kerr et all 
(120121) : the other five are given in Tables [4] and 0 which also 
list available flux density measurements. 

PSR J1747-4036 was observed mainly at the NRAO Green 
Bank Telescope (GBT), using the GUPPI snectrometciFl to 
sample a bandwidth of 800 MHz centered at 2 GHz, with typ¬ 
ical integration times of 5 min. 

The remaining MSPs were observed exclusively at Parkes 
where we first used the analog filterbank/PMDAQ data acqui¬ 
sition system as employed in the discovery observations, and 
later a digital filterbank (PDFB3/4), in all cases centered at 
1.4 GHz. Each observation typically lasted for 1 hr, except 
for the brighter PSRs J1902-5105 and J1903-7051 (about 
20 min). Apart from a dense set of observations to obtain or¬ 
bital parameters for the binary pulsars, and to unambiguously 
establish pulse numbering, we aimed to detect each pulsar ap¬ 
proximately monthly over a span of 1-2 yr. Especially for 

15 https://safe.nrao.edu/wiki/bin/view/CICADA/GUPPiUsersGuide. 
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Figure 2. PSR J1514—4946 and its variable flux density due to scintillation in the ISM, observed at Parkes at 1.4 GHz with 256 MHz of bandwidth using the 
analog filterbank system. All profiles are shown with 64 phase bins per period (P = 3.6 ms), with a linear greyscale, and with arbitrary phases between plots, 
(a) The best detection among 79 observations spanning 2 years, showing the signal strength gently increasing during the observation, with the integrated profile 
at the top. (b) An extremely weak detection, obtained only by folding 6 hr of data modulo the period predicted from the timing solution, (c) Signal strength as 
a function of radio frequency for the observation also shown in (a), (d) A 0.7 hr good detection with average flux density half of that shown in (a) and (c), and 
with very different time-integrated frequency structure, (e, f) Two observations separated by 9 hr (of length 1.0 hr and 1.2 hr, respectively), showing very different 
frequency structure. In many instances, one or two features of ~ 10 MHz in width dominate the detected signal. Often, no such feature is present and the pulsar 
is extremely faint, or not detectable. 


PSR J1514-4946, this required a very large number of ob¬ 
servations, because the pulsar is very faint on average and its 
received radio flux density varies enormously due to interstel¬ 
lar scintillation (Scction |2.2. 1 1 ). We improved the precision of 
the timing solutions by substantially extending the measure¬ 
ment baselines with a few additional observations in 2014 and 
early 2015. 

We used the PRESTO and PSRCHIVE dHotan et al.ll2004l) 
software to analyze the raw data and obtain pulse times-of- 
arrival (TOAs). For the MSPs with sparse radio detections, 
we obtained initial estimates of the orbital parameters using 
the method described in iFreire et all ( 2001 ). We then used 
TEMPC[3 and TEMP02 dHobbs etakl 2006 ) to obtain phase- 
connected timing solutions. Starting with solutions spanning 

1 yr for each pulsar, we detected gamma-ray pulsations for 

16 http://tempo.sourceforge.net. 


six MSPs. For PSR J1903-7051 we then obtained a few LAT 
TOAs. While the radio TOAs have much higher precision, 
they span only 3 yr. Adding gamma-ray TOAs increases the 
solution span and improves some resul ts, in particular the 
proper motion measurement (see Section [2.5.2l) . 

2.4. Polarimetry 

In principle, study of the polarized radio emission from pul¬ 
sars can constrain the magnetic field geometry and, particu¬ 
larly when considered together with gamma-ray profile char¬ 
acteristics, can elucidate emission locations and mechanisms 
(cf. Section [T2l) . 

In order to measure the polarization characteristics of each 
pulsar, we used the digital filterbank PDFB3 at Parkes, and 
GUPPI at GBT in coherent dedispersion mode, to do a few ob¬ 
servations recording calibrated folded full-Stokes pulse pro¬ 
files. These data were analyzed in standard fashion with 
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Table 3 

Preliminary Parameters for Four Millisecond Pulsars 



PSR J0955-6150 

PSR J1012-4235 

PSR J1036-8317 

PSR J1946-5403 

Right ascension 3 , R.A. (J2000.0).. 

09 h 55 m 39 s 

10 h 12 m 07 s 

10 h 36 m 20 s 

19 h 46 m 24 s 

Declination 3 , deck (J2000.0). 

-61° 48'36" 

-42°35'01" 

-83°17'01" 

-54°02'46" 

Spin period, P (ms). 

1.999 

3.101 

3.408 

2.710 

Dispersion measure, DM (pc cm -3 ) 

160.7 

71.6 

27.0 

23.7 

Orbital period, (d). 

24.578 

37.972 

0.335 

0.130 

Projected semi-major axis, x (1-s) . 

13.283 

21.263 

0.506 

0.0435 

Eccentricity, e . 

0.11 

<0.001 

<0.001 

< 0.001 

Companion mass' 3 , m 2 (Mq). 

>0.21 

>0.26 

>0.14 

>0.021 

Galactic longitude, / (deg). 

283.7 

274.2 

298.9 

343.9 

Galactic latitude, b (deg). 

-5.7 

11.2 

-21.5 

-29.6 

DM-derived distance 0 , d (kpc).... 

3.8 

2.5 

1.0 

0.9 


Note. — The listed P and DM values are from the discovery observations. Orbital parameters are from fits to sets of Doppler-shifted barycentered periods. 
a These discovery pointing positions have ±7' uncertainties. On the assumption that the MSPs are associated with the target gamma-ray sources, in some cases 
the position is better constrained (see Table [2}. 

b Derived from the pulsar mass function f\ assuming m\ = 1.35 Mq IThorsett & Chakrabarty|[T999f) and i < 90°. f\ = /P{j) 2 Tq = (m 2 sin*) 3 /( m i +ra 2 ) 2 , 
where Tq = GM pJc 3 = 4.9 25 //s, m\ and m 2 are the pulsar and companion masses, respectively, and i is the orbital inclination angle. 
c Using the ICordes & Laziol l2f>02 l Galactic free electron density model. The individual estimates have substantial uncertainties. 


PSRCHIVE. The resulting profiles are shown in Figures0]and 
0] and rotation measures (RMs) are listed in Tables 0] and [5] 

The main profile component of PSR J1514-4946 is about 
80% linearly polarized, with a flat position angle (PA.) of 
linear polarization (Figure^). 

PSR J1658-5324 has a complex profile, with a high (up to 
100%) linear polarization fraction in most pulse components 
(Figure 0]t). Such a high level of linear polarization is rela¬ 
tively uncommon among MSPs (see lYan et al.l20fHlDai et al.l 
120151) . Despite the PA. being measured over a large span of 
rotation phase, we were not able to obtain a sa tisfactory rotat¬ 
ing vector model fit (iRadhakrishnan & Cooke! 19691) . 

The profile of PSR J1747-4036 is curious. This MSP has 
the second largest value of DM/P among those known in the 
Galactic disk, and so was particularly well-suited for observ¬ 
ing with a coherent dedispersion system. We show two such 
observations in Figure^- At 1.5 GHz, the profile has two fea¬ 
tureless components, with a very flat PA. throughout. Fack of 
variation in PA. across the entire profile is unusual (there is no 
such example in a well-studied sample of 24 MSPs; lDai et al.l 
l2015i ). The profile observed at 2 GHz (which overlaps in fre¬ 
quency with that at 1.5 GHz) looks similar. At 0.8 GHz, how¬ 
ever (bottom panels of Figure 0};), the second component is 
brighter than the first, signifying that it has a steeper spec¬ 
trum, and a third component appears (more easily discernible 
in linear polarization), with P.A.s offset by about 90° from the 
rest of the profile. At this frequency, emission is detectable 
across 85% of the pulse period (but the observed profile may 
be somewhat affected by multi-path propagation; Cordes & 
Fazio 2002 predict pulse broadening of 0.02 P). 

The PSR J1902-5105 profile is relatively unusual in show¬ 
ing no discernible linear polarization (there is possibly a very 
small amount of circular polarization in both principal com¬ 
ponents; Figure 0}l). 

The radio profile of PSR J1903-7051 is shown in the bot¬ 
tom panel of Figure 0] It displays a weak, highly-polarized 
component (pulse phase (f> ss 0.58) leading a much stronger 
double-peaked component (0.65 < (j) < 0.85). The first peak 
is highly linearly polarized and also shows a degree of circular 
polarization, while the trailing (brightest) peak is completely 
unpolarized. 


2.5. Gamma-ray Observations and Analysis 

The gamma-ray properties of the first five MSPs we dis¬ 
covered (PSR J0101-6422 and those in Table 01 ) have been 
reported elsewhere dKerr et al.l 1201 2t lAbdo et afll2013l) . and 
we focus here on the newly-detected gamma-ray pulsations 
of PSR J1903-7051. To characterize its properties, we ana¬ 
lyzed “reprocessed Pass 7” SOURCE class Fermi -LAT events 
collected between 2008 August 4 (the start of nominal oper¬ 
ations) and 2015 February 1, and with energies between 0.1 
and 30 GeV. The data were filtered to exclude events whose 
reconstructed direction exceeds a zenith angle of 100° and 
those taken when the observatory was rocked more than 52° 
from the zenith, or when the region of interest (ROI) around 
the pulsar approached the Earth’s limb. These cuts minimize 
the bright gamma-ray background contribution from the limb 
of the Earth. 

2.5.1. Spectral Analysis and Photon Weights 

We considered a 20°x20° region around PSR J1903-7051, 
and used the binned likelihood formulation of gt- 
like (Fermi Science Tools v. 09-35-0203) to fit an 
exponentially cut-off power law to the pulsar emis¬ 
sion, dN/dE = No(E/Eo)~ r exp(-E/E c ). We used the 
P 7REP_S0URCE_V15 model of the instrument response 
function^ and the isotropic and diffuse background models 
of the 3FGF catalog (lAcero et all [2015b . as well as the 
point sources therein. The best-fit parameters and their error 
estimates are listed in the last section of Table 0] With models 
for the pulsar and background sources, we can compute 
a weight w, for each photon i giving the probability that 
the photon originated from the pulsar (iKeril 12011 al) . These 
weights allow improved separation of the pulsar signal from 
its background in the following analyses. 

2.5.2. Pulsar Timing 

Folding FAT data for PSR J1903-7051 based on the radio- 
only timing solution results in a clear drift in the position 
of the gamma-ray profile peak outside the time interval cov¬ 
ered by the radio ephemeris. Thus, using the s ame se lec- 
tion criterion as for the pulsar light curve (Section 12.5.31 ), we 

17 http://fernii.gsfc.nasa.gov/ssc/data/analysis/software. 

18 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html. 
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Table 4 

Parameters for Four Millisecond Pulsars with Coherent Timing Solutions 



PSR J1514-4946 

PSR J1658-5324 

PSRJ1747-4036 

PSR J1902-5105 


Timing parameters 

Right ascension, R.A. (J2000.0) . 

15 h 14 m 19? 1141(1) 

16 h 58 m 39?34359(9) 

17 h 47 m 48?71692(3) 

19 h 02 m 02? 84821(9) 

Declination, decl. (J2000.0). 

-49° 46' 15 7 516(5) 

-53°24'07"003(1) 

-40°36'54"773(1) 

-51°05'56"9695(8) 

Proper motion in R.A., dcos(<5) (masyr -1 ). 

-0.3(32) 

0 . 2 ( 8 ) 

- 0 . 8 ( 6 ) 

-4.8(13) 

Proper motion in decl., S (mas yr -1 ). 

-30.0(66) 

4.9(23) 

-4.9(16) 

-4.4(16) 

Spin frequency, / (Hz). 

278.60300920296(5) 

409.95436264371(4) 

607.67753906573(2) 

573.92104496683(5) 

Frequency derivative, / (Hz s -1 ). 

-1.4473(8) x 10 - 15 

-1.8746(6) x 10 - 15 

-4.8510(5) x 10 - 15 

-3.0301(4) x 10 ^ 15 

Epoch (MJD). 

55520.0 

55520.0 

55520.0 

55520.0 

Dispersion measure, DM (pc cm -3 ) . 

31.05(2) 

30.81(3) 

152.98(1) 

36.25(1) 

Orbital period, P/ y (d) . 

1.922653523(5) 



2.0118037388(9) 

Projected semi-major axis, x (1-s) . 

1.933268(2) 



1.9019570(7) 

Time of ascending node, T asc (MJD) . 

55585.8605555(3) 



55162.2815604(1) 

esino; a , EPS1 . 

6.453587(2) x 10 ~ 6 



5.5239429(7) x lO ^ 6 

ecos oj a , EPS2 . 

8.789469(3) x 10 ~ 6 



-1.9671264(9) x 10 - 6 

Span of timing data (MJD) . 

55160-57013 

55166-57057 

55161-56993 

55161-57014 

rms timing residual (^s) . 

8.3 

2.9 

2.4 

3.8 

Flux densities ’ 3 and rotation measures 

0.8 GHz flux density, So 8 (mly) . 



4.8 


1.4 GHz flux density, S| 4 (mJy) . 

0.17 ±0.11 (TV = 71) 

0.7 ±0.5 (N = 37) 


1.2 ±0.5 (TV = 69) 

1.5 GHz flux density. Si 5 (mJy). 



0.9 


2 GHz flux density, S 2 (mJy). 



0.5 ±0.1 (TV = 28) 


Rotation measure, RM (rad m -2 ). 

35 ±15 

4 ± 7 

-39 ±2 


Derived parameters 0 

Spin period, P (ms). 

3.589 

2.439 

1.645 

1.742 

Characteristic age, r c (10 9 yr). 

4.9 

3.5 

2.0 

3.1 

Spin-down luminosity, E (10 34 erg s _1 ). 

1.0 

3.0 

11.3 

6.7 

Surface dipole magnetic field strength, B (10 8 G) 

2.1 

1.7 

1.5 

1.3 

Eccentricity 3 , e . 

( 1.1 ± 0 . 2 ) x 10~ 5 



(5.9 ±0.7) x lO ^ 6 

Mass function, f\ (Mq). 

0.00210 



0.00183 

Companion mass, m 2 (Mq). 

>0.17 



>0.16 

Spectral index 4 , a . 



R -2.5 


Galactic longitude, / (deg). 

325.25 

334.87 

350.21 

345.65 

Galactic latitude, b (deg). 

6.81 

-6.63 

-6.41 

-22.38 

DM-derived distance, d (kpc). 

0.9 

0.9 

3.4 

1.2 

Composite proper motion, // (mas yr -1 ) . 

30.0 ± 6.6 

4.9 ±2.2 

5.0± 1.6 

6.5 ±1.4 

Transverse velocity, Vj_ (km s _1 ) . 

R 130 

fa 20 

R 80 

R 40 


Note. — All uncertainties are reported at the 1 a level. Numbers in parentheses represent the TEMP02 timing uncertainties on the last digits quoted. 
a The orbital eccentricities were derived from the parameters (e sincere cos lj), fitted using the TEMP02 ELL1 binary model (Lange et all2Q0lt) . 
b Si .4 and S 2 values are averages and standard deviations for N detections. Individual flux densities were determined by computing the area under each pulse 
profile compared to its off-pulse rms, scaled using the measured system equivalent flux density at the location of the pulsar. Values for PSR J1747-4036 are 
from GBT observations; So 8 and Si 5 are from single flux-calibrated observations (Figure [3j) and have uncertainties < 10%. For one similar observation at 
2 GHz, S 2 = 0.45 mJy. 

c The following have been used: r c = P/(2P), E = An 2 x 10 45 P/P 3 erg s _1 , and B = 3.2 x 10 19 (PP) 1 / 2 G, with P in s, where P = 1//. The listed values of these 
parameters include corrections for acceleration effects (mainly due to proper motion; see lCamilo et all 1994) . 
d S v oc v a , where S v is the flux density at frequency v. 


have extracted 16 gamma-ray TOAs following the method of 
iRav et all (1201 111 and used these to extend the timing solu¬ 
tion. The resulting parameters, including an improved mea¬ 
surement of the proper motion, appear in Table @ 

2.5.3. Light Cun’e 

The gam ma-ray light curve in the top pa nel of Figure 0] 
with H-test (Ide Jager et al.lll989tlKerril201 lbl) value of 231, is 
a probability-weighted histogram of the photon phases. Be¬ 
cause the probability weights w,- allow an optimal “soft” cut, 
no tuning of the photon selection is required, and the his¬ 
togram includes all events within 2° of the pulsar position. 
The structure of the light curve is robust against additional 
background, so here, unlike in the spectral analysis, we do 
not apply the cut on zenith angle when the ROI is near the 
horizon, increasing the total number of photons by ~ 25% 


and slightly increasing the background level. Error bars for 
a bin follow the typical prescription a 2 = w h with the 
sum over the N ly photons in the / 1 h bin. The background 

level is given by (XS vv ' _ zC; w ?)/-^bin, with the sums over 
all Nry photons in the profile and ;V| lln the total number of 
bins. This level, based on the spectral model, represents the 
expected contribution from all diffuse and background point 
sources. Both pulsed and unpulsed emission from the posi¬ 
tion of PSR J1903-7051 will show up as a signal in excess of 
this background level. However, from inspection of the light 
curve, it is evident that PSR J1903-7051 has no substantial 
unpulsed component. A slight excess may be due to an in¬ 
creased background over the spectral model due to the less 
stringent zenith cut. 
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Figure 3. Polarimetric pulse profiles for PSRs J1514—4946 (a), J1658—5324 (b), J1747-4036 (c), and J1902—5105 (d). In the bottom sub-plots, the black line 
corresponds to total intensity, while the red and blue traces represent linear and circular polarization, respectively. In the top sub-plots, the position angle of linear 
polarization (P.A.) is plotted for bins in which the linear signal-to-noise ratio > 3, corrected to the reference frame of the pulsar using the RMs listed in Table[4] 
PSR J1747-4036 was observed at the GBT with GUPPI (recording a bandwidth of 800 MHz at a central frequency of 1500 MHz and 200 MHz of bandwidth at 
820 MHz, with the two profiles aligned by eye). All others were observed at Parkes with PDFB3 recording 256 MHz of bandwidth. 


2.6. X-ray Observations 

In an attempt to assist in determining accurate positions 
for the pulsars before timing solutions existed, we under¬ 
took X-ray observations of the fields of PSRs J1514-4946 
and J1658-5324 with the ACIS-S camera on board the 
Chandra X-ray Observatory ( CXO ). For both observations 
the (then) best pulsar positions were centered on the back- 
illuminated S3 chip, and we analyzed the event data with 
the s tandard CIAO version 4.7 software (Fruscione et al. 
120061) . After excluding events outside the energy range 
0.3-7.0keV we searched the fields for any point-like X- 


ray source possibly associated with the pulsars using the 
cell detect source search tool. In each case we detected 
a faint source 0E3 ± 0."6 away from the pulsar timing posi¬ 
tion. We name these, respectively, CXOU J151419-494615 
and CXOU J165839-532406, which we identify as the pulsar 
counterparts on the basis of positional coincidence. Observa¬ 
tion and source parameters are listed in Table [6] 

For three other MSPs, Swift X-Ray Telescope (XRT) obser¬ 
vations have been undertaken as part of a Swift campaign of 
observations of Fermi -LAT unassociated sources. No sources 
were detected at the pulsar positions. More information, in- 
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Table 5 

Radio and Gamma-ray Parameters of PSR J1903—7051 


Parameter Value 


Timing parameters 

Right ascension, R.A. (J2000.0). 19 h 03 m 38?7935(3) 

Declination, decl. (J2000.0). -70°51'43"461(2) 

Proper motion in R.A., dcos((5) (mas yr -1 ). -8.8(16) 

Proper motion in decl., 8 (mas yr -1 ). -16(2) 

Spin frequency, / (Hz). 277.94006243351(8) 

Frequency derivative, / (Hzs -1 ). -8.06(4) x 10 -16 

Epoch (MJD). 56526.0 

Dispersion measure, DM (pc cm -3 ). 19.66(1) 

Orbital period, P b (d). 11.05079833(2) 

Projected semi-major axis, x (1-s). 9.938869(2) 

Time of ascending node, r aS c (MJD). 56027.2292914(7) 

esincu, EPS1. 2.0261968(5) x 10" 6 

ecosu;, EPS2 . 1.1855635(5) x 1(T 7 

Span of timing data * * 3 (MJD). 54760-57013 

ims timing residual 3 (/xs). 12.8 


Flux densities 15 and rotation measure 

1.4GHz flux density, Si .4 (mJy). « 0.6 (.N = 9) 

3.1 GHz flux density, S 3.1 (mJy). ~ 0.14 (N = 3) 

Rotation measure, RM (rad m -2 ). 11 =L 24 


Derived parameters 0 

Spin period, P (ms). 3.597 

Characteristic age, r c (10 9 yr). 7.1 

Spin-down luminosity, E (10 33 erg s -1 ). 6.8 

Surface dipole magnetic field strength, B (10 8 G) 1.7 

Eccentricity, e . (2.0 =b 0.5) x 10 -6 

Mass function, f\ (Mq). 0.00863 

Companion mass, m 2 (Mq) . > 0.28 

Spectral index, a. . ~ —1.8 

Galactic longitude, / (deg). 324.39 

Galactic latitude, b (deg). -26.51 

DM-derived distance, d (kpc). 0.8 

Composite proper motion, p (masyr -1 ). 18.3 d= 2.0 

Transverse velocity, V± (km s _1 ). « 70 


Gamma-ray parameters 

Gamma-ray-radio profile lag d , 8 (P) . 0.55 d= 0.05 

Gamma-ray (> 0.1 GeV) photon index, T . 1.90 =b 0.16 

Gamma-ray cut-off energy, E c (GeV). 7.8=L4.0 

Photon flux (> 0.1 GeV) (10“ 8 cnr 2 s" 1 ). 1.2=L0.3 

Energy flux (> 0.1 GeV) (10 -11 erg cm -2 s -1 )_ 0.9 =b 0.2 


Note. — All uncertainties are reported at the 1 <7 level. Numbers in 

parentheses represent the TEMP02 timing uncertainties on the last digits 
quoted. Uncertainties for gamma-ray spe ctral parameters are statistical only 
(for a discussion of systematic errors, see lAcero et alJ2015h . 

3 We have used both radio and gamma-ray TOAs to derive this timing solu¬ 
tion. 

b The listed values are median flux densities for N calibrated detections. In¬ 
dividual S L4 values ranged over 0.13-1.5 mJy. In four 3.1 GHz observations 
we did not detect the pulsar. 

c The listed values of r c , E, and B include corrections for acceleration ef¬ 
fects. 

d This is measured from the profiles in Figure[4] with the gamma-ray profile 
centroid at phase 0 = 0.27 and the radio profile reference phase mid-way 
between its full observed span (0.55 < (j) < 0.9). 

eluding flux and luminosity upper limits, is given in Table [6] 
All these X-ray detections and upper limits are consistent 
with the known distribution of MSP X-ray luminosities (e.g., 
iPossenti et al.ll2002l and references therein). 

3. DISCUSSION 

In our directed radio survey with the Parkes telescope of 56 
unidentified Fenni-LXI gamma-ray sources, we detected 11 



0.0 0.2 0.4 0.6 0.8 1.0 

Pulse Phase 


Figure 4. Gamma-ray and polarimetric radio profiles of PSR J1903—7051. 
Top: Weighted >0.1 GeV light curve (6.5 yr of LAT data), displayed with 30 
phase b ins, inc luding estimates of the error bars and the background level (see 
Section 1753) . Bottom: Phase-aligned Parkes 1.4 GHz pulse profile (black: 
total intensity; red: linear polarization; blue: circular polarization). The green 
error bars are centered on the position angles of linear polarization (displayed 
for RM = 0 with an offset of +130° for ease of view). 

MSPs, 10 of them discoveries (Section [2Tt ). 

3.1. Ten New Millisecond Pulsars 

Two of the 10 new pulsars are isolated (PSRs J1658-5324 
and 11747-4036; Table [4j. Another five are in circular orbits 
(e < 10 -3 ) with 0.2-0.3 Mq companions, likely helium-core 
white dwarfs (PSRs J0101-6422, J1012-4235,11514-4946, 
J1902-5105, and J1903-705 l: lKerr et alJlMl and Tables^ 
Q] and |5}. Three of these have orbital periods of about 2 d, 
while the other two have periods of 11 d and 38 d. These char¬ 
acteristics are typical of MSPs known in the Galactic disk. 
The timing precision of PSR J1747-4036 is good enough that 
it has been added to the NANOGrav gravitational wave pulsar 
timing array (lMcLaughlinll2013l) . 

The remaining three systems (see Tabled are less common. 
PSR J0955-6150 is in a P h = 24 d orbit with a m 2 « 0.25 M 0 
companion, but with a very significant eccentricity (e = 0.11). 
It joins four Galactic MSP systems with broadly similar pa¬ 
rameters: 0.2 < m 2 < 0.3 M W , 0.03 < e < 0 .13, 22 < Pt < 
32 d (see lAntoniadisl20T4tlKnispel et all2015L and references 
therein). Such eccentric systems are not predicted through the 
stand ard MSP formation channels (see Phinnev & Kulkarni 
1994ft, leading to alternative scenarios dFreire & Taurisll2014t 
Anton iadis 20 01). 
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Table 6 

X-ray Observations of Five Millisecond Pulsars 



PSR J1514-4946 

PSRJ1658-5324 

PSRJ1747-4036 

PSR J1902-5105 

PSR J1903-7051 

Telescope/Instrument. 

CTO/ACIS-S 

CXO/ACIS-S 

Swij f/XRT PC 

Swifi/XRT PC 

Swifi/XRT PC 

Exposure (ks). 

9.9 

9.9 

3.3 

4.2 

3.2 

Background-subtracted counts. 

9 

23 

< 6.6 a 

< 12.7 a 

<6.6 a 

A H b (10 21 cm- 2 ). 

1 

1 

5 

1 

0.6 

X-ray flux c , /o.i- 2 . 4 keV (10 14 ergcirT 2 s _1 ) 

1.1 

2.9 

<23 

< 12 d 

< 7 

X-ray luminosity e , Lq. 1 - 2.4 keV (lO -4 #).... 

1.1 

0.9 

<28 

< 3 

< 8 


Note. — A blackbody model with kT = 0.2keV is assumed in all cases, absorbed by the indicated column density Ah. Considering instead a pow er-law 
spectrum with photon index T = 2, and calculating flux/luminosities for the 2-10 keV range, does not fundamentally alter our conclusions (Section l2.61 . 
a We used 47”-radius extraction regions around each pulsar (90% PSF radius for XRT), and 141” radii to estimate backgrounds. For two sou rces we obtained 
zero background-subtracted counts, while for PSR J1902—5105 we obtained three counts. In each instance we convert to a 3 c upper limit (see lGehrelslI 19861) . 
b Absorbing columns are estimated from the DMs according to /Vj|(l() 2l cm 2 ) = 0.03DM fsee IHe et ali2013h . 

0 We used http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html to obtain these unabsorbed flux estimates (or 3 a limits). 
d lTakahashi et ah (2013 ) report a somewhat lower flux limit based on a 38 ks Suzaku XIS observation. 

e Isotropic 0.1-2.4 keV luminosities are given in terms of the acceleration-corrected values of E, and for the DM-derived distances, listed in Tables[4]and[5] 


PSR J1036—8317 is in an 8 hr orbit with a ss O.16M 0 com¬ 
panion, similar to so-called “redback” systems (where out¬ 
flows from >0.15 M 0 non-degenerate companions can cause 
irregular radio eclipses of the pulsar; e.g., [D’Amico et all 
120011) . Despite many observations at essentially all orbital 
phases, no radio eclipses have been detected. This is there¬ 
fore unlikely to be a redback. More likely the companion is 
a white dwarf. If its distance is « 1.0 kpc, as inferred from 
the DM, this could be a good target for optical studies; it 
may be a system similar to PSR J1012+5307, an MSP in a 
14 hr o rbit with a spectroscopical ly identified 0.16 M 0 white 
dwarf (Ivan Kerkwiik et alJ fl996h - or to PSR J1738+0333 
(lAntoniadis et alJl2012l). 

PSR J1946-5403 is in a 3 hr orbit with a > 0.021 M 0 com¬ 
panion. These parameters suggest a “black widow” inter¬ 
acting system (sub-day binaries with degenerate < 0.05 M 0 
companions; e.g.. iFruchter et allll988l) . Most such systems 
show radio eclipses near superior conjunction, but so far we 
have not detected any (when folded using the known orbital 
parameters, the pulsar is detected in all six data sets listed in 
Table [T] including two observations of superior conjunction). 
This could be due simply to geometry: if the actual compan¬ 
ion mass is slightly larger than the minimum value inferred 
from the mass function, e.g., if mi > 0.025 M 0 , the orbital 
inclination angle i < 60° and we could be viewing the system 
relatively face-on. In any case, with a DM-derived distance of 
0.9 kpc, this may also be an interesting optical target. 

It seems curious that among the 10 MSPs discovered in 
this Parkes survey only one is in an interacting binary sys¬ 
tem (either black widow or redback), when about 50% of the 
67 MSPs so far discovered in searches of unidentified LAT 
sources worldwide are in such systems (see [Roberts! 12013 ). 
For instance, in a recent survey at Arecibo, five of six MSPs 
discovered are either black widows or redbacks (H. T. Cro- 
martie et al. 2015, in preparation). However, the Arecibo 
searches have integration times of 7 =15 min, which are far 
more suitable for the discovery of few-hour binaries than the 
1 hr Parkes integrations: the maximum accelerations probed 
by our searches scale as T~ 2 . At the GBT, where half of these 
67 MSPs were discovered, typical integration times are SO¬ 
TS min (e.g.. iRansom et all 1201 ll P. Bangale et al. 2015, in 
preparation). In any case, the pulsars that we did not detect 
unbiasedly in every one of our Parkes survey observations are 
the two sub-day binaries (T able [Qi, at least in part due to large 
and rapidly changing orbital acceleration (cf. FigureQ}- 


3.2. Six New Gamma-ray Millisecond Pulsars 

Radio timing observations since 2009 have yielded rota¬ 
tional ephemerides for 43 of the 67 MSPs discovered in LAT- 
guided searches. Of these, 39 are now confirmed as gamma- 
ray MSPs_(the other four are unrelated to the LAT sources; 
e.g.. iKeith et all 120111) . including the six discovered in our 
Parkes survey for which we already have timing solutions. In 
this paper we have for the first time presented the ephemerides 
and polarimetry for five of these MSPs, and gamma-ray prop¬ 
erties for PSR J1903-7051. In some cases these results add 
to our understanding of the pulsars summarized in the second 
LAT catalog of gamma-ray pulsars (2PC: lAbdo et al.l[2013l) . 

For PSR J1903-7051, the radio peak leads the gamma-ray 
profile by 6 ~ 0.5 P. Other pulsars with such values of <5 have 
gamma-ray profiles composed of only one peak (2PC), which 
is also the case for PSR J1903-7051 (Figure [4]). This is one 
of the few MSPs for which spectral curvature (deviation from 
a power-law spectrum) is not apparent even with 4 yr of LAT 
data (Table |2). Our spectral fits to 6.5 yr of data yield the 
largest cut-off energy of any pulsar, even if with large uncer¬ 
tainty: E c = (7.8 ±4.0)GeV (Table[5j see also FigurelSf). The 
MSPs with the next largest values of E c are also from our 
Parkes sample: PSRs J1747-4036 and 11514-4946 (2PC). 

With our proper motion measurement for PSR 11514-4946 
(Table [4ji, its intrinsic E is 1.6 times smaller than previously 
thought. Its computed gamma-ray efficiency, already high in 
2PC (30%), now increases to 50%. As usual this assumes a 
geometry-dependent beaming correction of /o = 1 (see 2PC 
for definition), i.e., assuming isotropic emission. Perhaps for 
PSR J1514-4946, /n<Cl. Or maybe its DM-derived distance 
of 0.9 kpc is an overestimate (its nominal transverse velocity 
is the largest of those listed in Tables [4] and 0 although it is 
not unusually large for an MSP; e.g. JGonzalez et akll201 ll) . 

PSR J1747-4036 is notable in several respects. Spectral 
curvature is not apparent in the 3FGL catalog (see Table [2]). 
Its E is large (fourth highest among the 2PC MSP sample). 
Its nominal distance (Table [4]) is also among the largest in 
the 2PC MSP sample, and it is a relatively faint gamma-ray 
source (10er in 3FGL; Table 0). If it had a typical E at that 
distance presumably it would not be a detectable LAT pulsar. 
This large E ultimately is due to its very short spin period 
(third shortest among disk MSPs). While PSR 11747-4036 
may be somewhat extreme in this regard, it does appear that 
LAT-detected MSPs are a smaller-P, higher-/-’ population than 
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radio-selected MSPs (iRav et alJl2012l) . In our Parkes sam¬ 
ple, PSRs J1902-5105 and J0955-6150 (the latter apparently 
quite distant) are also short-period MSPs, with P <2 ms. 

PSR J1902-5105 is one of only six MSPs known with 
phase-aligned gamma-ray and radio light curves (2PC). This 
points to emission that is either extended, and caustic in na¬ 
ture, or originates near the neutron star surface dVenter et alj 
12012 ). In order to constrai n their emission and viewing geom¬ 
etry, Johnson et al.l (120141) have jointly modeled the gamma- 
ray and radio profiles of 40 LAT-detected MSPs. They 
fit gamma-ray light curves using standard outer magneto¬ 
sphere gap (OG) and two-pole caustic (TPC) geometric mod¬ 
els assuming a vacuum retarded d ipole m agnetic held. For 
PSR J1902-5105, [Johnson et alj (120141) consider altitude- 
limited (al) versions of the standard models as well as a low- 
altitude slot gap model (laSG). They find that all three models 
match the gamma-ray and radio profiles fairly well. How¬ 
ever, the complete lack of polarization that we observe (Fig¬ 
ure [3]J) strongly favors the alTPC and alOG models, where 
high-altitude c austics are predicted to have a strong depolariz¬ 
ing effect dDvks et alj|2004l) . over the laSG model, where one 
expects non-zero polarization from near-surface emission. 

For PSR J1514-4946, the high level of observed linea r po- 
larization and modest PA. swing (Figure[3ji and Section [ZAt 
suggest we are viewing the edge of a cone beam at relatively 
low altitude. For PSR J1658-5324, the very high level of lin¬ 
ear polarization suggests that at least some components orig¬ 
inate at relatively low altitude (the small peak with low po¬ 
larization, which might also originate from the opposite pole, 
could be a high-altitude caustic; Figure^). PSR J1747-4036 
has a very unusual polarization pattern (Figure [3};), which 
does not fit either standard radio core/cone or high-altitude 
caustic emission models. For these MSPs, iJohnson et all 
(120141) find models that plausibly match the gamma-ray light 
curves, but the radio profiles are not well reproduced. 

3.3. Sun’ey Statistics 

Our survey of unidentified LAT sources had a success 
rate of 20% (11 MSP s detected unbiasedly in 56 sources 
searched; Section 12.11 ), despite significant selection effects 
(Section [T2l i. This is very encouraging — a substantial num¬ 
ber of unidentified LAT sources contain previously unknown 
MSPs that can be detected at Parkes, and our source selec¬ 
tion criteria have allowed us to target them with satisfactory 
efficiency. H oweve r, comparing the initial segment ofthe sur¬ 
vey (Section [2TH) with the latter portion (Section [2.1.21 ) re¬ 
veals a drop in efficiency: we had to search 3.5 times as many 
sources in 2012, for twice as long in the aggregate, to discover 
as many pulsars as earlier on. 

Since our observations, another six pulsars have been dis¬ 
covered in these very same sources: four via direct pulsation 
searches of the gamma-ray photons, one of which was subse¬ 
quently detected as a very faint radio source with the GBT (the 
others might be significantly affected by scintillation, or per¬ 
haps they may yet prove to be gamma-ray-only pulsars); one 
at the Giant Metrewave Radio Telescope (GMRT) at a lower 
frequency than at Parkes (perhaps the Parkes non-detections 
reflect a steep radio spectrum); and one that is a transient ra¬ 
dio source — not active when we did our Parkes searches, but 
now detectable there (see notes l-o in Tabled. Considering 
these additional detections, a full 10 of the original 14 targets 
of our survey are now known to contain pulsars! This reflects 
superb source selection criteria for those targets. Why has 
the success rate decreased since then? The original (1FGL- 


based) sources were brighter in gamma rays, allowing for rel¬ 
atively unambiguous determination of spectral characteristics, 
on which we based our target selection. There might be sec¬ 
ondary effects related to this: e.g., the brighter lFGL-based 
targets might be nearer to the Earth on average, and any as¬ 
sociated radio pulsars could then also be brighter on average, 
although we do not see such an effect among our small sample 
of detected MSPs. 

3.4. Spectral Characteristics of LAT Pulsars 

Now that the third catalog of LAT sources is in hand, based 
on substantially more data and a much better understanding 
of the instrument and background than was available ear¬ 
lier, we can usefully revisit our 56 sources and the MSPs 
found amid them. In Table [2] we have summarized some 
properties of these sources as obtained from the 3FGL cat¬ 
alog. Given that gamma-ray pulsars typically have exponen¬ 
tially cut-off power-law spectra (lAbdo et al.l 12013l) . by con¬ 
trast to power-law spectra for AGN, it is not surprising that 
the spectral curvature gleaned from 3FGL correlates well with 
gamma-ray pulsars: of the 12 pulsars with known gamma- 
ray pulsations listed in Table [2] (the eight classified as “PSR,” 
as well as J1035.7-6720, J1227.9-4854, J1624.2-4041, and 
J1744.1-7619), 10 display significant curvature according to 
3FGL (Curve > 4cr; TableO. However, two other established 
gamma-ray MSPs (PSRs J1747-4036 and J1903-7051, re¬ 
ported here) show no curvature, according to 3FGL (Curve < 
3 a). Now that 170 gamma-ray pulsars are known, it is to be 
expected that some may depart from the norm. In addition, 
“Curve” in 3FGL tests against a log-parabolic spectral model, 
which may not be a good proxy for testing against exponen¬ 
tially cut-off power laws, especially when the source is not 
very bright or the background is problematic (as it happens, 
these two MSPs are the faintest in gamma rays of the 12 un¬ 
der consideration; see “Sig” in Table [21)• Thus, in considering 
whether a particular 3FGL source is a good pulsar candidate, 
we should do more than simply look for significant cataloged 
curvature. 

3.5. Spectral Classification of Unidentified LAT Sources 

Flux information is available from the 3FGL pipeline pro¬ 
cessing in five energy bins (two within 0.1-1 GeV, two within 
1-10 GeV, and one for 10-100 GeV). After familiarizing our¬ 
selves with how known gamma-ray pulsars (and control 
sources) appear in five-bin spectra by visual inspection of 
many 3FGL plots, we have classified the 56 target sources 
of our Parkes survey according to a heuristic ranking scheme 
that we describe next. 

We qualitatively assess the likelihood of a source being a 
pulsar. According to this scheme, a classification of “1” de¬ 
notes near certainty of being a pulsar; “2” is less conclusive 
but quite plausible; “3” is a poorer pulsar candidate; “4” is 
very likely not a pulsar; “5” is not a pulsar. 

These classifications are based on combinations of the fol¬ 
lowing spectral characteristics (both are listed for each source 
in TableO: strong energy cutoff (c); parabolic (peak at center 
energies, p); flat or rising spectrum at < 1 GeV (r) — c, p, r 
are “positive” features leading to higher likelihood of a source 
being a pulsar (for examples, see Figure[5j: and e). Power law 
(1); variability (v); monotonically decreasing (d) or increasing 
(i); excess high energy emission (flat or rising spectrum; h) — 
1, v, d, i, h are “negative” features, not ordinarily associated 
with pulsars (see Figure[5ji and d). Capitalization means that 
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Figure 5. Examples of spectral classification of unide ntified Fermi-L XT sources with a view towards identifying pulsar candidates. The panels show the spectral 
energy distribution (SED) from the 3FGL analysis JAcero et all2015b , with the dashed lines/gray bands indicating the best-fit model SED (either a power law or 
a curved log-parabolic spectrum) and its uncertainties. Discrepancies between the spectral measurements a nd th e model indicate that the fit may be unreliable. 
Our classifications (1-5) and qualitative spectral characteristics that underlie them are described in Section 13751 The full list of 56 classified sources searched 
in our Parkes radio survey is given in Table[2] (a) An example of a poor pulsar candidate, consistent with a power law (1), with both high-energy emission (h) 
and an AGN association; (b) an example of an atypical pulsar spectrum, now associated with the state-changing MSP J1227-4853; (c) an excellent candidate 
with a strong cutoff (C), a spectrum slightly more peaked than that of the typical pulsar (P), and rising at < 1 GeV (r); (d) a poor candidate spectrum with both 
high-energy emission (h) and unreliable spectral points (?); (e) a good candidate with an unusually peaked spectrum (P); (f) a decent candidate with evidence for 
a cutoff (c) but an unusually steep low-energy spectrum, monotonically decreasing (d); now known to be PSR J1903—7051 (this work; e.g., see Tableland top 
panel of Figure[4). 


the features are more certain, except for p/P which indicates 
the sharpness of the parabola. A “?” indicates poor spectral 
quality, increasing uncertainty in classification. An asterisk 
represents an odd spectrum with “banana” shape to high en¬ 
ergy (see Figure [5j>). In Figure 0 we show 3FGL spectral 
plots for six of our Parkes targets that illustrate the features 
described above. 

Sixteen of our 56 sources are classified as 5 or 4 (or 3 with 
a blazar association), and we regard them as no longer viable 
pulsar candidates. 

Of the 12 known gamma-ray pulsars in Table [2] nine are 
classified as 1 and two others (already noted in Section [3741 
as not curved in 3FGL) as 2. The one confirmed gamma- 
ray MSP classified as a 3, PSR J1227-4853, is a spectral 
outlier (Figure [5J)). In a rare state-changing binary system, 
it is borderline variable within 3FGL, and recently displayed 
significant variability (see iJohnson et al.ll2015l). This and its 
sister system J1023+0038 (IStap pers et all 20141). as well as 
the young pulsar .12021+4026 (lAllafort et al .1120131) . are coun¬ 
terexamples to the usual assumption that pulsars are steady 
gamma-ray emitters. Of the other five MSPs known in Ta¬ 
ble [2] (none of which has yet a reported long-term timing so¬ 
lution or detected gamma-ray pulsations), one is classified as 
1, one as 2, and three as 3. Based on prior statistics, we ex¬ 
pect that most, perhaps all, of these five MSPs will eventually 
be established as gamma-ray pulsars. That most have a rel¬ 
atively poor classification in our scheme likely reflects their 
faintness (three, all discovered in our survey, have 3FGL sig¬ 
nificance of 6.6-8.8 a), or location close to the Galactic plane 


(PSR J1536-4948 is at b = 4.8°). 

3.6. Further Searches of LAT Sources 

Following from the discussion in Section l331 while we cer¬ 
tainly recommend additional searches first of the “1” sources 
that remain without coincident pulsars (of which there are 
seven in Table O, followed by the “2” sources (nine in Ta¬ 
ble [2j, the “3” sources (seven in Table O are also reasonable 
candidates for further searches^ Among these 23 sources, 
two are coincident with globular clusters that currently have 
no known pulsars, and a population of MSPs in those clusters 
is a plausible origin for the gamma-ray emission. We expect 
further MSP discoveries among these 23 targets: while the 
largest offset between our optimized search and the 3FGL po¬ 
sitions is 3', which is not very large compared to the 7(2 half¬ 
width at half-maximum Parkes beam, many promising targets 
have not been search ed e nough to counter the selection effects 
presented in Section |T21 — e.g., four of the seven “1” sources 
have been searched at Parkes only once or twice at locations 
not far from their 3FGL positions (see Tablc[2]i. In addition, in 
the near future we intend to re-analyze our existing data sets 
with “jerk searches,” i.e., accounting for changing accelera¬ 
tions that are relevant for short orbital periods (cf. Figure|T}. 

After all such searches are done, it remains a possibility that 
through a combination of faintness and extreme orbital and 
spin parameters, some of those promising LAT sources could 

15 One of the “2” sources, 3FGL J1417.5-4402, has recen tly been reported 
to be a binary system currently with an accretion disk IStrader et al.ll2015l ). 
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still harbor undetected radio MSPs beamed towards the Earth. 
It is also possible that some of these sources may be MSPs that 
are detectable only via their gamma-ray emission, although 
this fraction is known to be small (e.g jRomanill2012h . 
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